OPTIMIZED SPACE SHUTTLE TRAJECTORY SIMULATION 
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TRAJECTORY STRUCTURE 
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*The state equations (to be discussed later) define the position, velocity, and weight of the vehicle. 


STATE EQUATIONS 
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general 


For some problems it is necessary to be able to simulate a trajectory which branches at one or more points. For example, 
the optimization of a space shuttle ascent trajectory in which constraints are imposed during the booster return flight 
and during possible abort return flights (left half of figure) requires the capability to simulate a trajectory with multiple 
branching points. After each ascent trajectory simulation, the normal booster return flight and hypothetical abort trajecto- 
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TYPICAL ASCENT TRAJECTORY WITH TRAJECTORY SECTIONING AND GENERAL 

BOOSTER ABORT/RETURN CONSTRAINTS BRANCHING CAPABILITY 




GENERAL MATHEMATICAL PROBLEM 
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PERFORMANCE INDEX 







«TtAI & 

+ 

XI 

tH ^ 

■T— J 
•<— ( 

X 



X) 






0 






pP 












fcxD 

•2 






T3 






0 

cl 






P 

•rH 






4J 






PH 






O 






CO 






0 






p 



u 



r*H 



CD 



P 



H-J 



> 


4J* 

0 

CJ 



r-H 


rP 

s 



P 


bO 

P 



g 


•i— i 

CD 

Pi 

P 



■rH 

<4-H 


£ 





H-> 

CD 

bD 

<-< 



S 


P 

.9 



r— j 



H- 3 

CO 


P 


e» 

U$ 


•pH 


• 

CD 

s 


H- 3 
• pH 


bJD 

rP 

Pc 

CO 

.§ 


• 

CD 

g 

#4 

• 

fan 

• 

0 

•4-H 

0 

o 

s 

F-i 

CD 

H-> 

CD 


§3 

JD 

§ 

0 

Iti 


-j- 3 

g 

s 


fl 


CD 

O 

Pc 

0 

8 

P 

2 

.2 


p 

•+J 

tH 

P< 



ro 

0 

0 

Pc 

Pi 

P 


s 

0 

a 

o 

PC 

J“ H 

P 


Pc 

Pi 

Pc 

> 


O 

P 

C 

0 

0 

0 


*4-C 

•pH 

Oh 


fi 

P 

£ 

0 


CD 

Q« 

Pf 

O 

0 

rH 

O 

• pH 
rH 

H- 3 

2 



c -4— 1 

s 

2 

p 


O 

•rH 

Pi 

H- 3 

Ph 

0 

P, 

ad 

Pc 

P 

p 

ft 

0 

T3 


CD 

r"j 

o 

a 

X 

c-i 


S 

•pH 

<-H 

tw 

4-H 

•- 


rf 

2 

O 

o 

T3 


C<j 

a 

P 

1 

Pc 

0 

pQ 

g 

CD 

’g 

§ 

*o 




P 

2 


e 

c 

«*< 

P 

c 

CO 

o 

CD 

0 

0 

CD 

-p» 

o 

•pH 

H- 3 

s 

£ 

6 

3 

.& 

Ph 

O 

0 

CQ 

w 

CO 

CO 

o 

CO 

•pH 

•pH 

•pH 

•pH 

CO 

he 

P 





Pc 

*, E 


-©- 

i-3 

0 

Ph 

•S 

T5 

1-3 

P 

rt 

P 

P 





CO 

O 

CD 

F-< 




0 

6 

Ph 

CO 

0 

5 . 

CD 




T3 

rH 

Pc 

-g 




2 

O 

£ 




p 

a 


1445 



EQUATIONS FOR CONSTRAINT CATEGORIES 
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Integral inequality (e.g. specification of a maximum total heating parameter): 
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NONLINEAR PROGRAMMING TECHNIQUE 
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none of the inequality constraints are violated. A penalty term is added to the performance index, which keeps the 
solution in the feasible region at each stage in the iteration. 
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Each iteration requires minimization of the function P for a given value of r. This problem is termed the P-problem, and 
its solution is fundamental to the method of Fiacco-McCormick. The computational algorithm which was implemented 
minimizes the function P, using a Fletcher-Powell^ function minimization algorithm in conjunction with the simultaneous 
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GENERATION OF A FEASIBLE SOLUTION - This section reviews the algorithm developed by A. V. Fiacco for gener- 
ating initial feasible solutions. An arbitrary initial design vector y° will result in general in s< m of the inequality con- 
straints being satisfied while the remainder will be in violation. With no loss in generality, assume the constraints to be 
reordered so that the first s constraints are the satisfied constraints. Constraint s + 1 can be brought into the feasible set 
while simultaneously guaranteeing that constraints 1, 2,...,s remain in the feasible set by solving the following nonlinear 
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CONSTRAINED OPTIMAL CONTROL AS A NONLINEAR PROGRAMMING PROBLEM 
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PSEUDO-STATE VARIABLES 
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GENERATING A FEASIBLE SOLUTION FOR STATE AND CONTROL VARIABLE 
VARIABLE INEQUALITY CONSTRAINTS 12 
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however, that the integration limits can be updated at the end of each one dimensional search. Such frozen limits on the 
first term are not required. 



A DYNAMIC INEQUALITY CONSTRAINT TRAJECTORY 
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VIOLATED REGION 



TRAJECTORY SIMULATION 
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DIFFERENTIAL EQUATIONS OF MOTION 



1459 


= the horizontal component of the gravitational acceleration vector (positive in 
local due north direction) 







a .a 

£ a 

>* 5 

co O 
<D co 


Cd 


03 

G 

g ^ 
C 2 
o ^ 
o .5? 
° 13 

>» » 
4-» <D 

•3 * 

r\ Cd 


id 

> 


•G 03 

^3 M 

(D 


ID 

H-S 
M o 

d ^ 

G <d 

£ o 

£13 


03 

.S 

O 

o 

o 


<d 

o 

03 

a 

00 

o 

«H 

<d 

03 


id 

*G 

4 -» 4 -» 

"G vh 

o o 

'o G 

> O 

• i— i 

id +5 

> GO 

'■£ o 

rS ^ 
'o +2 
»-H G 

_ CD 

* b 

■S 3 

T 3 2 
(D *G 
X ^ 

■c +3 

o ro 


•a w> 

.2 ‘g 

L- W? 


bD 

G 

■^H 

ID 

T3 


CD 

^3 

(D 

*o 


Cd 


% S 

■§ <4 § 

4-* r- 1 , 03 

'§!•§ 

Ilf 

>, E 

ti u i_ 

'o >; S 

O +j 3 

> 2 £ 
£ <L> ^ 

.5 xl 

•ti 1 3 c 

2 -a «* 

g c .a 

H Cd +J 
» ~ >. 
5 ^ al 

e j{ 4g 

c3 U> . 

c > 

VO "* •£ 

*3 6 ^ 

§ o,ii 


in 


t* <D 
1 ZJ <u 

fj »! 
i ■£ -S 

S cd ^ 
i 2 13 « 
id ^ G 


(D 

G 


Q 

ID 

> 

CD 

4— 1 

03 


ID 

G 


rs 1-1 

k 0 

- 4-4 o 

o 5 

ID > 

> 

>* +-» 
+2 cd 

4-» 
co 


O > 

CO <L) 

' ' > 


o 

5 <> 

Si x! 


G id 

CO -4 

5 ,SP 


S O 

2 c 
o 0 
x; *-1 a 

® £ I 


g 
< u 

c — 

g “ 

13 QJ 


>> % 2 
4-4 <D G 

O X +- 

° 1 o 

E<^« 
1 3 

CD co 

jG > 4-» 


^ +D Qji 
+-H _, 

o £ 


ID 

4-4 

cd 

o g C 

O o £ 

CO »i— j ^ 

CO 4-4 Q_> 

03 cd co 


CD 


CD 


I 


CD 


G 

s -se 

O 


<D 

<^6 


CD 

G 


CD >> 

*3) .2 

§0 

| > O 
2 « «r 

^ > <u 

13 S 

CS 


•a 

c „ ^ 

« 

^ Ji 

o X J3 

nl +-• 
+j > T3 

at > G 

t, 3 

og2 

<D .SP 

n * O 

cd 4-4 
G cd 
G G 

h S 2 
• S 

» 8 

S g 

CD G 

b cr 

^ .1 8 
v *- 1 CTl 

M si 

Cd C*_i 
G O ID 
W CD .5 

i — | Tl I»h 

S ' a •§ 
.. S .2 

s -g» « 

“ 0 ^ 3 
U l— 1 

S^s 

«.3 * 

-5 ID 

a g’S 

e 1 2 

cd 2 cd 
G > 

CD ~ 

CD 


‘O. • 


CD 

1 ■ t 

W) 


ID 

G 

G^ «-» 

^ <cr G 

. ID 

rs 

X)<w ® 

s^- 

<u <U 
£ X « 

t*H 

U -S 

■o iJ o 
o o> 
co -g D< 

G co 

CD CD 

^ ^ 

T3 

^0 


G<^ 

4—* 


G - 

o O 

cd «-< v 


cd 


G 

o 

cd 


(D o3 

’3) ^ 

G O 
& ^ 


caX) 

C 

cd 


G 

o 


G Oi 
O 0 
G 


G 

C 

cd 

G 

ID 

G 

4-4 

G 

G 

O 


G 

bD 

G 

O 

t-4 

G 


b£) 

G 

cd 

cd 

<D 

G 

4-4 

G 

bO 

G 

O 

•— 

G 


G <£T 

4-» W' 

.25 co 

g x 

cc3 cd 

45 -g 


X 

cd 


O 

t-H 

CD 

G 

4-4 

G 

G 

G 

O 

«-i 

cd 

G 

O 


CD 

G 


G 

O 

G 

cd 

G 

O 


« .a 


G 

ccS G G 

cd 

JJ . 

O 

G G „ 

(D O > 

cd<-H 1 . . 

^_» ’W' C+_( 

O 


T3 

G 

O 

O 

CD 


O 

1-1 


cd 

G 


cd 


x q 

cd .G 


G 

bO 


CD 

G 


.S3 ^ 


>, 0 


X 

cS 

C3 

<D 

Xi 


3 

O 

X> 

at 

G 

O 


X! 

(30 

‘C 

tD 

X 


G 

O 

X 


G 

O 


cd 

+-» 

O 


<3 < < 'E 

>. 

■ 4-4 

^ (N cn < 


C3 


o 

8 

b<J 

.K 


"O 


53 


5> 


K 

8 


CD 

■« 


■8 


1460 


L_ 



COORDINATE SYSTEMS AND ATTITUDE ANGLES 




The central body total gravitational acceleration vector (g) is obtained directly from the central body gravitational poten- 
tial energy function (U) by taking its gradient (VU). The expression for U is shown in Equation 25. 

The assumed gravity model includes up through the fourth harmonic term (N = 4) and assumes symmetry about the 
rotational axis. The corresponding expression for the gravitational components is g r = 9U/9r and = -(9U/90)/r. 
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propellant flow rate are constant to models in which the vacuum thrust, propellant flow rate, and engine throttling coeffi- 
cient are tabular functions of generalized independent variables. Of particular use to the space shuttle ascent simulation is 
a model which automatically provides throttling to meet thrust acceleration constraints. For all models, the effective 
thrust is computed as a function of the atmospheric pressure. 




SPACE SHUTTLE BOOSTER RETURN EXAMPLE 
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Angle of attack > 0 deg. 

Angle of attack s 90 deg. 

The last two constraints were due to the unavailability of aerodynamic data outside of this region. 



AND C, AS FUNCTIONS OF MACH NUMBER AND ANGLE OF ATTACK 
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The initial angle of attack and bank angle were each modeled by 12 parameters as a function 
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